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Light hadron spectroscopy 

D. K. Sinclair^*t 

'^High Energy Physics Division, Argonne National Laboratory, 
9700 South Cass Avenue, Argonne, Ilhnois, 60439, USA 

I review the progress that has been made in light hadron spectroscopy from lattice QCD, since the LATTICE'94 
conference in Bielefeld. 
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1. INTRODUCTION 

Calculation of the hadron spectrum has been a 
central goal of Lattice QCD for 20 years. A preci- 
sion calculation of the hadron spectrum would be 
a validation of QCD, and a validation of Lattice 
QCD as a calculational technique. It would give 
us confidence when we apply lattice QCD to cal- 
culating other quantities. In addition we would 
be able to clarify the nature of observed hadrons 
and predict where others might be found. True 
precision hadron mass calculations have been per- 
formed for heavy-quarkonium systems and 
are beginning to appear for heavy-light systems 
g. Despite some monumental efforts such as that 
by the GF-11 group ||], the hght hadron spec- 
trum remains considerably less well determined 
than one might hope for. This is unfortunate, 
since this part of the hadron spectrum is the most 
complex and least well understood. 

The reasons for this situation are well known. 
The light hadrons (with possible exception of the 
glueballs) are large, their size being determined 
by the pion Compton wavelength, so that large 
lattices are needed. Small lattice spacings are 
needed to make contact with the continuum limit. 
The u and d quark masses are small, making the 
Dirac operator ill-conditioned and expensive to 
invert. There are a large number of excited states, 
which makes it difficult to isolate the ground state 
contribution to (almost) any hadron propagator. 
At small quark masses, hadron propagators be- 
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come more noisy - especially for large (time) sep- 
arations. Simulations with dynamical quarks suf- 
fer from critical slowing down as the quark mass 
approaches zero. 

For these reasons, much of the work on light 
hadron spectroscopy has been performed in the 
quenched approximation. Unfortunately, it is in 
just this regime of light quarks that this approx- 
imation becomes the most suspect. Even here, 
few measurements have used lattices of much 
more than 3.5fm in spatial extent. Many calcu- 
lations are restricted to u and d quark masses 
so large that m^/mp > 0.5, requiring consid- 
erable extrapolation to obtain physical results 
(mTr/wp ~ 0.18). 

Contributions since LATTICE'94 are listed in 
Tables 1-8. In each table, the columns give p, 
lattice size, number of configurations and number 
of source positions/configuration. 

In addition, three groups, COLUMBIA, S. Kim 
- DKS and Gottlieb have contributions on the 
chiral limit. Kilcup has presented calculations of 
the connected and disconnected contributions to 
the rj' propagators using staggered quarks. 

Before proceeding to summarize the above con- 
tributions, let me first make a few general obser- 
vations on directions pursued since LATTICE'94. 
For work prior to this one should consult recent 
review articles such as those by Michael |^ and 
by Weingarten 0, and the collections of parallel 
talks (and poster sessions) in the proceedings of 
LATTICE'94 and earlier conferences in the series. 

In the past year the trend has been towards 
high statistics, rather than larger lattices and/or 
smaller quark masses. Most of the work has been 
in the quenched approximation. Many of these 
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GERMAN: 


6.0 16^ 


X 32 
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LANL: 


6.0 32^ 


X 64 


150 




JLQCD: 


6.0 24^ 


X 64 


200 
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JLQCD: 


6.0 24^ 
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5.7 16^ 
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Table 4 

QUENCHED VALENCE 
Liu - Dong: 6.0 163 x 24 



Table 5 

IMPROVED-QUENCHED CLOVERLEAF 



SCRI: 



6.80 
7.40 
7.60 
7.75 
7.90 



163 X 32 
83 X 16 
163 X 32 
163 X 32 
83 X 32 



Table 6 

IMPROVED-QUENCHED IMPROVED 

CORNELL: 6.8 53 x 14 
7.1 63 X 16 



Table 7 

FULL QCD STAGGERED 

COLUMBIA: 5.7 163 x 40 



Table 3 

QUENCHED CLOVERLEAF 



UKQCD: 



5.7 123 X 24 480 
6.0 163 X 48 125 
6.2 243 X 48 184 



Table 8 

FULL-QCD-STAGGERED CLOVERLEAF 
SCRI: 5.6 163 x 32 100 
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light hadron spectrum calculations were per- 
formed as an integral part of other projects, such 
as matrix-element calculations and/or heavy- 
light physics calculations. 

This high statistics, sometimes combined with 
use of multiple sources/sinks to isolate the ground 
state, has enabled the spectra to be measured 
with sufficient precision to permit extrapolations 
in quark mass and lattice spacing. Tests have 
been made as to how well one reproduces the 
spectrum of hadrons containing strange valence 
quarks. 

Using Wilson quarks for spectrum calculations 
has the advantage that flavour symmetry is ex- 
act, although this is at the expense of losing all 
chiral symmetry. A perhaps more serious limita- 
tion is that Wilson quarks introduce errors 0{a). 
(a is the lattice spacing.) Staggered quarks have 
the advantage of keeping a vestige of chiral sym- 
metry, at the expense of losing flavour symmetry. 
While having the disadvantage that their lattice 
spacing is effectively 2a they have the advantage 
that their errors are 0{a^). An alternative solu- 
tion is to use an improved Wilson-quark action, 
the Sheikholeslami-Wohlert |^ or cloverleaf ac- 
tion, which has the flavour symmetry of the Wil- 
son action but errors ©(a^). 

Preliminary studies aimed at implementing the 
Lepage programme |^ of using improved actions 
to allow use of larger lattice spacings and smaller 
(in number of sites) lattices have been under- 
taken. 

Section 2 presents highlights of the above men- 
tioned spectrum calculations. In section 3 we 
discuss the chiral limit of the pion spectrum of 
quenched lattice QCD. Section 4 gives a brief dis- 
cussion of the work by Kilcup et al. on the rj' 
mass, while section 5 is devoted to discussions 
and conclusions. 

2. HIGHLIGHTS OF SPECTRUM CAL- 
CULATIONS 

What follows are brief descriptions of re- 
cent contributions on light-hadron spectroscopy, 
which have been summarised in Tables 1-8. Most 
of the figures in this section were provided by the 
people mentioned, or prepared by me from data 



provided by said persons. All work presented in 
this section should be considered as preliminary. 
The light-hadron spectroscopy contributions of 
two of the JLQCD projects was only a sideline 
to other calculations, and was not available to 
me at the time of preparing this paper. I there- 
fore refer you to the parallel talks by Hashimoto 
0] and Aoki ^iov details (I thank M. Okawa 
for bringing these works to my attention). 

2.1. QUENCHED STAGGERED 

Gottlieb has calculated the light hadron 
spectrum with staggered quarks in the quenched 
approximation to lattice QCD. He used "corner" 
wall sources and point sinks. His lattice sizes 
ranged from 8^ x 48 to 32^ x 64. At /? = 5.7 
and P — 5.85 his quark masses ranged from 0.01 
to 0.16, while at /3 = 6.15 they ranged from 0.005 
to 0.08. 

High statistics on multiple lattice sizes has 
enabled him to analyse finite size effects. For 
m^r/mp > 0.5, finite size effects appear to be 
fairly small for spatial boxes of extent > 1.5-2fm. 
Use of several masses at each beta has enabled 
him to extrapolate mtqjmp to wiq = (a linear 
extrapolation appears adequate). Finally, use of 
several /3's and hence a's has permitted extrapo- 
lation to a = 0. Both a linear and the accepted 
quadratic extrapolations in a give acceptable fits 
and mass ratios which are too high by 5-8%, 
which is normally considered acceptable consid- 
ering the systematic uncertainties of the ratios, 
fits and error estimates, (see Figure |l|). 

A plot of mN /nip versus rriT^/mp for these sim- 
ulations is shown in Figure ^. 

Gottlieb ||l^ has also attempted to calculate 
the quantity J, introduced by Lacock and Michael 
p^ ] as a consistency check on the meson spec- 
trum. (J — mic evaluated at mv/mps = 

mji* /rriK-) His very preliminary results are: 
0.56(1), for P = 5.7 (24^ x 48 lattice), 0.55(2) for 
P = 5.85 (20^ X 48 lattice) and 0.49(1) (24^ x 48 
lattice), and 0.45(1) for /3 = 6.15 (32^ x 64 lat- 
tice). A quadratic extrapolation to a = gives 
0.43(1) compared with an experimental value of 
0.48-0.50. He has also studied the PCAC rela- 
tionship for the pion mass (see Section 3). 

S. Kim and D. K. Sinclair ||l^ have studied the 
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Figure 1. a dependence of m-^ /vrip^ showing linear 
and quadratic fits. 



Figure 2. Edinburgh plot - quenched staggered. 



light hadron spectrum with quenched staggered 
quarks at /3 = 6.0 (their previous simulations had 
been at /3 = 6.5 10). They used wall sources and 
point sinks. Lattice sizes of 16"^ x 64, 24'^ x 64, 
and 32'^ x 64 were used to study finite size effects. 
Quark masses of 0.0025, 0.005, 0.01 were used. 

High statistics enabled us to study finite size 
effects. The tt showed no significant finite size ef- 
fects for box sizes > 24^^; similarly for the p (at 
least for 2 heavier quark masses). Because of the 
difficulty in finding a plateau in the nucleon ef- 
fective mass on the 32^ x 64 lattice, and because 
of the persistent difference between the two dif- 
ferent estimates (sources) for the nucleon mass, it 
remains unclear how significant finite size effects 
are for the nucleon for these box sizes and quark 
masses. 



As one can see from the effective mass plots 
for the nucleon (Figure ^ , at least for the largest 
lattice, it is unclear whether a plateau is attained 
before the signal is swamped by noise, despite the 
high statistics. We thus concluded that a wall 
source was a poor choice for a 32"^ x Nt lattice at 
[i = 6.0. 

From the extrapolated p mass, a^^ = 
1.87(2)GeV at /? = 6.0. This yields ae.o/ae.s = 
2.02(3) compared with scaling prediction {MS) 
of 1.8976. Finally we include an Edinburgh plot 
(Figure The 2 points for each mass/lattice 
size come from the 2 different nucleon sources. A 
rough estimate of the J parameter, is J w 5.6, 
compared with the experimental value of 4.8-5.0. 
Our results for the chiral behaviour of the theory 
are presented in section 3. 

S. Kim and S. Ohta are performing 
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Figure 3. Nucleon (A^i) effective mass plot - quenched staggered. 



quenched simulations at /3 = 6.5 on a 48^ x 64 lat- 
tice. On their current sample of 50 lattices they 
have calculated the light hadron spectrum with 
staggered quark masses 0.01, 0.005, 0.0025 and 
0.00125 (in lattice units) enabling them to probe 
down to m^/mp < 0.3. They use wall sources 
and point sinks. 

When combined with the earlier work of Kim 
and Sinclair [|l6|, this work should enable them 
to study finite size effects at /3 = 6.5, where scal- 
ing is apparent and violations of flavour symme- 
try (indicated by the difference in masses between 
Goldstone and non-Goldstone pions) are minimal. 
Their preliminary results, summarized in the Ed- 
inburgh plot of Figure || would seem to indicate 
considerable finite size effects in the rriN /mp mass 
ratio. A word of caution is appropriate, since 
these are relatively low statistics, preliminary re- 
sults. mN/mp decreases both because tun de- 
creases and because rup increases. Since the low 



statistics forces them to fit propagators at rela- 
tively small time separations, at least some of the 
difference could be due to not having yet reached 
the plateau in their propagators. Only higher 
statistics can resolve this issue. 

2.2. QUENCHED WILSON 

The German collaboration (G. Schierholz et 
al.) have performed quenched light hadron spec- 
troscopy with Wilson fermions as part of a nu- 
cleon structure function calculation [Q. They 
used a 16'^ x 32 lattice at /3 — 6.0, and used a 
Jacobi smeared source and sink. 

Using high statistics with Jacobi smearing 
(which appears to be a good choice), they ob- 
tained well behaved hadron propagators with 
good evidence for plateaux (see Figure ^ which 
shows the nucleon effective mass). They also 
use 2-component "non-relativistic" sources (pro- 
jected with i(l + 74)), which halves the number 
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Figure 4. Edinburgh plot - quenched staggered. 



Figure 5. Edinburgh plot - quenched staggered. 



of required inversions. The hadron mass values 
they obtain are given in Table 9. 

The LANL collaboration (communicated by 
R. Gupta) have performed hadron spectroscopy 
with quenched Wilson quarks on a 32'^ x 64 lat- 
tice at /5 6 [0. They use four k's — 0.153, 
0.155, 0.1558, 0.1563 — for u, d, s physics (they 
also do some charm physics). Hadron propaga- 
tors are obtained for wall source/point sink (WL), 
Wuppertal source/point sink (SL) and Wupper- 
tal source/ Wuppertal sink (SS). Since SL and SS 
effective masses approach the infinite separation 
limit from above, WL from below, a linear combi- 
nation was used to reduce the contributions from 
excited states. 

No significant finite size effects were observed 
between 24^ x Nt and 32^ x 64 lattices. Lin- 
ear extrapolation to the correct pion mass gave 
mN/mp = 1.42(4), niA/mp = 1.80(5), (both 10 
- 15 % too high) and a"! = 2.319(44)GeV (from 
nip). Figure |^ shows the linear dependence of 
meson and baryon masses on the quark mass. 



Figure || is the Edinburgh plot from these sim- 
ulations. 

The values of required to give the correct 
values for mx/mp^ uiK'/TTrip, and m^/mp differ 
by ^ 20%. Even with the largest of these es- 
timates for m^, the baryon octet and decuplet 
splittings are still as much as ~ 30% smaller than 
experiment. 

Several forms for the lattice dispersion relation 
were tested for mesons. The best was found to 
be sinh^(i;/2) = sin^(p/2) sinh^(m/2), which 
worked well up to charmonium masses. 

The JLQCD collaboration (presented by 
A. Ukawa) have simulated quenched lattice QCD 
on a 24^ x 64 lattice at /3 = 6.0 H]. They have 
calculated hadron propagators and hence masses 
for 3 different k values of Wilson quarks. Most 
of their preliminary results were for k — 0.1550 
where m-n/mp « 0.70. 

They have used their high statistics to exam- 
ine the "wiggles" in the plateaux of effective mass 
plots. Such behaviour is common. It is worrisome 
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Figure 6. Nucleon effective mass v's temporal 
separation. 



because these fluctuations are often considerably 
larger than 1 standard deviation, and thus unex- 
pected. 

What Ukawa and his colleagues do is to cal- 
culate the complete covariance matrix from the 
data. They then model the data with a prop- 
agator which is a 2-exponential fit to the data, 
with gaussian fluctuations in the directions of the 
eigenvectors of the true covariance matrix. The 
width of these gaussian fluctuations is determined 
by the eigenvalues of this true covariance matrix. 
Thus this model has the same average propaga- 
tor and the same covariance matrix as the actual 
data. Its effective mass plot for any finite sample 
shows deviations from the plateau by more than 
1 standard deviation, as does the real data. They 
thus conclude that the existence of such large fluc- 
tuations is a property of the covariance matrix 
of the theory, rather than an indication that no 
plateau exists or that their datasets are not inde- 
pendent. 



Table 9 

The hadron masses in lattice units at /3 = 6.0. 

K 

0.1515 0.153 0.155 

0.504(1) 0.422(2) 0.297(2) 

nip 0.570(2) 0.507(2) 0.422(2) 

TUN 0.900(4) 0.798(4) 0.658(4) 



Suzuki and collaborators |2l| have studied 
quenched QCD with Wilson quarks on a 16^ x 32 
lattice at /3 = 5.7. Their study is aimed at deter- 
mining whether abclian monopolcs, which have 
been suggested by 't Hooft j22| to be responsible 
for confinement, are responsible for generation of 
hadron masses in the chiral limit. 

Their procedure is as follows. After gauge fix- 
ing to the maximal abelian gauge, they factor- 
ize the gauge fields into a diagonal abelian gauge 
field, and a residual field. They then calculate 
field strengths corresponding to these diagonal 
gauge fields. These abelian field strengths are 
separated into "smooth", "photonic" contribu- 
tions and Dirac string, monopole contributions. 
Abelian gauge fields are reconstructed from each 
of these field strengths. 

They have calculated Wilson quark and hence 
hadron propagators for each of the full gauge 
field, the abelian gauge field, the monopole (Dirac 
string) gauge field and the photonic gauge field. 
What they have found is that the tt and p masses 
for the full and abelian gauge configurations ap- 
proach the same values in the chiral limit. The 
results for the proton are unclear. More work 
needs to be done before it will be known whether 
the masses for the monopole gauge field yield the 
same results. The "photon" field case, as ex- 
pected, does not reproduce the correct spectrum 
in the chiral limit. 

2.3. QUENCHED CLOVERLEAF 

I will start by discussing briefiy the work of La- 
cock and Michael jl^, who introduced the quan- 
tity J — niK' evaluated at my / raps = 
TTLK'/raK ~ 1.8. Experimentally J=0.49(2), if 
they use the assumption that m?pg oc (TOi+m2)/2 
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Figure 7. Dependence of hadron masses on the 
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Figure 8. Edinburgh plot - quenched Wilson, 
LANL and GERMAN collaborations 



with TOi and m2 the quark and antiquark masses 
in the pseudoscalar meson, PS and my = c + 
d{mi + 7712) with c and d constants. Ignoring 
for the moment any isospin breaking, m-jr/'Tnp 
fixes the degenerate u and d quark masses, while 
TUK/rnTr fixes the strange quark mass. Hence, on 
the lattice, J is determined. From the UKQCD 
data for the unimproved cloverleaf action they 
found J = 0.37(3) (4), and similar values for Wil- 
son data from other groups. 

The UKQCD collaboration (presented by 
H. Shanahan |Q) have studied the light-light 
hadron spectrum with tadpole-improved clover- 
leaf quarks. They ran at /3 = 5.7, 6.0 and 6.2 
on lattices sizes up to 2A^ x 48. Two k values 
were used at each (3 {k = 0.14144 and 0.14226 
at /3 = 6.2). Point and Jacobi smeared and/or 
"fuzzed" sources and sinks were used. 

Using multiple sources and sinks allowed sta- 
ble 2-particle fits to propagators and better de- 
termination of the ground state mass. They used 



smeared-smeared (SS), smeared-local (SL) and 
local-local (LL) propagators in their fits. Effec- 
tive mass plots for each of these propagators for 
the vector states (from their earlier runs with 
an unimproved cloverleaf action) are shown in 
Figure || jl^ . Note that these various propaga- 
tors clearly have different couplings to the excited 
states, which allows stable 2-mass fits (shown in 
the figure). 

J (defined above) was measured. The results 
were 0.368(5) (/3 = 5.7), 0.330(19)(17) (/3 = 6.0) 
and 0.363(60) (/3 = 6.2), still below the experi- 
mental value. 

my — « 0.55GeV^ experimentally, 

from light mesons through light-charm mesons. 
UKQCD 's measured values of this quantity are 
in good agreement with this observation. Wil- 
son quarks and unimproved cloverleaf quarks are 
in considerably worse agreement with this obser- 
vation. This is of particular interest, since this 
splitting is due to the spin-orbit coupling. 

Figure |l^ shows the a dependence of mp (ob- 
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Figure 9. Effective mass v's time separation for 
vector mesons, at k — 0.14266 The symbols are 
LL(+), SL(x) and SS(0). The curves are the 
2-state fits. 



tained by Hnear extrapolation in niq). What one 
notices is that the points tend to flatten out as 
one approaches a = 0, which is what is expected, 
since corrections are 0{a^) for the tadpole im- 
proved cloverleaf action. This contrasts with the 
linear dependence on a for the Wilson action. 

2.4. QUENCHED VALENCE 

K.-F. Liu and S.-J. Dong have studied light 
hadron spectroscopy in the valence quark approx- 
imation 1^^. They work on a 16^ x 24 lattice at 
(3 = 6, with K = 0.148, 0.152 and 0.154. 

A word of explanation is needed since the 
quenched approximation itself is often referred 
to as the valence quark approximation (even in 
this talk). Liu and Dong's valence quarks are 
quenched Wilson quarks with all backward mov- 
ing contributions removed. Hence for meson 
propagators, exactly 1 quark and 1 antiquark pass 
through any time slice between source and sink. 
For baryons exactly 3 quark lines pass through 



Figure 10. a dependence of p mass. K is the 
string tension. 



any timeslice. 

Light hadron spectroscopy with valence quarks 
shows even better SU(6) symmetry than one 
might expect. The A and nucleon masses are 
almost degenerate, as are the p and pion masses 
(in SU(6) A and N are in the same multiplet, as 
are p and tt). This is shown in Figure |ll|. 

They have also applied this approximation to 
the calculation of low energy matrix elements. 

2.5. IMPROVED-QUENCHED 
CLOVERLEAF 

The SCRI collaboration (presented by R. Ed- 
wards are calculating the light hadron spec- 
trum in the quenched approximation. Gauge con- 
figurations are being generated using the tadpole 
improved action of Lepage et al. [^0, with 
f3 — 6.8-7.9 (note that these /3 values correspond 
to much smaller /3's for the Wilson action). Quark 
propagators are being calculated using use tad- 
pole improved cloverleaf action. Lattice sizes are 
8^ X 16, 8^ X 32 and 16^ x 32. They are using 
two gaussian smeared sources (and sinks?). Four 
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Figure 11. A,N,p and tt masses as a function of 
quark mass in the valence quark approximation. 



k's corresponding to fixed pion masses are being 
used. Fits to multiple correlation functions al- 
low 2-exponential fits and better isolation of the 
ground state. 

This gauge action has errors 0{a*) (and 
©(a^a^)), compared with 0{a^) for the standard 
(Wilson) gauge action. The cloverleaf quark ac- 
tion has errors O(a^) (and 0{a^a)). They are 
measuring the a dependence of the spectrum. 
The meson spectrum appears consistent with be- 
ing linear in (the nucleon mass has yet to be 
measured). The a dependence of the p mass is 
given in Figure [l^. 

The lattice spacings are relatively large (0.16 - 
0.41fm). (The chief advantage of using the im- 
proved action is to allow this use of coarser lat- 
tices.) Of course, one really needs an O(a^) quark 
action to make full use of this fact. 

J = rriK' j"T , defined above, is measured and 



is < 0.39 for the parameters of these simulations 
(experiment = 0.48(2)). The potential and string 
tension are calculated as a scale for the hadron 




Figure 12. a dependence of the {p) mass, a is the 
string tension. 



2.6. IMPROVED-QUENCHED 
IMPROVED 

The Cornell group (presented by M. Alford 
p7[|) has studied quenched QCD with the tad- 
pole improved action of Lepage et al. . The 
quark and hence hadron propagators are gener- 
ated on these configurations using the D234 im- 
proved quark action. They work with 5'^ x 14 and 
6'^ X 16 lattices at /3 = 6.8 and 7.1 respectively. 

The improved gauge action produces results 
which differ from the continuum theory by terms 
O(a^) (and 0{ala^). The D234 quark action is 
given by 

6 

whose physics differs from the continuum by 
terms 0{a?') (and 0{asa?))- This allows the use 
of much larger lattice spacings than was possible 
with the Wilson gauge action, and Wilson, clover- 
leaf or staggered quark action. In fact, with the 
chosen /3's the sizes of the spatial boxes are ^ 2fm, 
which is larger than a 16^ spatial box at /3 = 6.0, 
with the Wilson gauge action. 

The spectrum appears similar to that obtained 
from conventional simulations - mt^jvnp ~1.5 - 
1.6 for mTj/mp '^0.6 - 0.65. The value of the 
quantity J, defined above, is 0.38(2) at /? = 6.8 
and 0.40(2) at /? = 7.1 
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2.7. FULL QCD STAGGERED 

The COLUMBIA group (presented by D. Chen 
|28| ) have continued their studies of the hght 
hadron spectrum with 2 flavours of dynamical 
staggered quarks. They ran on a 16"^ x 40 lattice, 
at /3 = 5.7, with staggered quark mass niq = 0.01. 
Two sources were employed for spectrum calcu- 
lations, a 16^ (wall) source and an 8*^ (octant) 
source, with the gauge field fixed to coulomb 
gauge. 

Using 2 sources, with high statistics, enables 2- 
particle fits for the pion propagator and 4-particle 
fits for other hadron propagators, enabling better 
isolation of the ground state in each channel, and 
giving a preliminary estimate for the mass of the 
first(?) excited state. 

For the A they tried Landau as well as 
Coulomb gauge. The mass from the Landau 
gauge propagators was smaller than that from the 
Coulomb gauge propagators. What is unclear is 
whether this is a true effect due to the unphysi- 
cal excitations in the Landau gauge, or just that 
the onset of the plateau is delayed in the Landau 
gauge. Figure |l^ shows the A effective mass in 
coulomb gauge. 



The chiral behaviour of the pion propagator as 
a function of valence quark mass has been calcu- 
lated and will be discussed further in section 3. 

2.8. FULL-QCD-STAGGERED 
CLOVERLEAF 

The SCRI group (presented by J. Sloan ^) 
are calculating the light hadron spectrum on the 
16^ X 32 HEMCGC configurations which were 
generated with dynamical staggered quarks. The 
spectrum is being calculated with tadpole im- 
proved clover leaf quarks. (3 = 5.6, and = 0.01 
for the dynamical quarks. Their spectrum calcu- 
lations use wall and shell-smeared sources. 

They employ fits to multiple propagators which 
enables better isolation of the ground state. 

Extrapolations of hadron masses to the chiral 
limit are being studied. Vector meson masses are 
fit to forms which are polynomial in the pseu- 
doscalar meson mass Mps- They have tried fits 
to polynomials up to fourth order, but which do 
not include a linear term. Preferred fits were to 
the forms 



(2) 



and 



My = Co + C2Mj>s + CiMj>s 



(3) 



" CO 

D O 

E 

a. ^ 

(D ^ 

t 6 



16''x40 /S = 5.7 ma = 0.01 1 src 



5 .10 15 

distance 
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The cubic fit above yields coefficients Cq = 
0.393(4) and C2 = 1.10(3) for cloverleaf quarks 
The data and these fits are given in Figure 14 
The fits for Wilson quarks are given for compar- 
ison. From these chiral extrapolations of the p 
mass they obtain = 2.25(3)GeV for Wilson 
quarks, = 1.92(2)GeV for cloverleaf. This is to 
be compared with 1.80(2)GeV obtained for stag- 
gered quarks. 

A linear approximation to J defined above, 
yields 0.425(6)(10) for Wilson quarks and 
0.444(5) (26) for cloverleaf compared with the ex- 
perimental value 0.499. Cloverleaf light vec- 
tor meson masses are in reasonable agreement 
with experimental values and closer than Wilson 
masses. Higher mass mesons are predicted to be 
lighter than experimentally measured. 



Figure 13. A effective mass in Coulomb gauge 
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aMy vs. aMpg^ 
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Figure 14. Cubic fit of my v's rripg 



3. THE CHIRAL LIMIT CHIRAL 
LOGARITHMS v's FINITE SIZE v's 
FINITE LATTICE SPACING 

In quenched QCD, cliiral perturbation theory 
predicts that 



rather than the usual 



(4) 



(5) 



||30| , pi|| where S is estimated to be « 0.2. 

On the other hand finite lattice size effects 
would be expected to modify the above relation 
(Equation 5) to 

,2 



= constant + 2iimq 



(6) 



where constant — > as the spatial volume V — > 
oo. The COLUMBIA group point out that these 
two departures from the canonical behaviour can 
be very difficult to distinguish . 

In our 32"^ x 64 pion mass data we observed 
departures from the normal relation above, which 



were well fit by the quenched chiral log formula 
1^ . Since there was good agreement between 
our 32'^ X 64 and 24'^ x 64 data, we rejected the 
interpretation of this as being a finite size effect. 
We show this data in Figure |l^. 
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Figure 15. nri^/niq v's at /3 = 6.0 (quenched). 



The COLUMBIA group first observed similar 
departures from the canonical relationship in full 
QCD at /? = 5.48 on a 16^ x 32 lattice and at 
(3 — 5.7 on a 16^ x 40 lattice when the dynam- 
ical quark mass was held fixed (at = 0.004 
and 0.01 respectively) Q. Although one would 
expect anomalous chiral behaviour in this case, 
it should not be identical to the quenched case. 
They also performed quenched simulations a,t (3 = 
5.7 on a 16^ x 40 lattice and obtained similar re- 
sults. For all their pion spectra they found good 
fits to the finite-size-modified formula. They also 
managed to fit our results to the formula with a 
finite intercept, although the fit had a much worse 

■ For their own simulations, they were able to 
correlate these finite size effects in to^ with those 
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in {ipip) . Figure |l^ shows the variation of jniq 
with quark mass. 




0.000 0.005 0.010 0.015 



Figure 16. m^/mq v's ruq. Columbia results 
(32^ X 64 results are from S. Kim and DKS). 



Gottlieb's quenched data at /? = 5.7, 5.85 and 
6.15, combined with a collage of data at /3 = 6.0, 
indicates rather different variation of /m^ with 
rUq at different /3's suggesting that at least some 
of the variation is a finite lattice spacing effect 
|l2| . His data at (3 = 5.7 shows mass dependence 
out to masses for which finite size effects must be 
small (see Figure |l^) . 

4. THE r;' MASS 

Kilcup et al. |Q have calculated the 77' mass 
using a method similar to that used by Kura- 
mashi et al. but using staggered, rather 

than Wilson quarks. They use 16^^ x 32 lattices 
at /3 = 6.0 (quenched) and /3 = 5.7 (full QCD - 
Columbia group). The full QCD configurations 
had mq = 0.025 and 0.01. 

They calculate both the connected (1-loop) 
and disconnected (2-loop) contributions to the 77' 
propagator. The connected piece is calculated in 
the standard fashion, but using a noisy source. 
The disconnected piece requires the inverse of the 
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Figure 17. m^/niq v's niq for quenched QCD - 
various /3's. (The (3 = 6.0 data at higher masses 
is essentially flat.) 



propagator at "zero" separation averaged over all 
sites of the source and sink time-slices (actually 
pairs of time slices with staggered quarks) . Since 
this is impractical, it is replaced by a stochastic 
estimator, using 96 noise vectors per lattice. 

The ratio of the 2-loop to 1-loop contributions 
to the propagator rises linearly in time separa- 
tion t for quenched QCD. For full QCD it ap- 
proaches a constant as t — s- 00, exponentially. 
Such behaviour is evident in their results. They 
estimate the rj' mass from the slope of the linear 
rise in the quenched theory and from the mass 
scale in the exponential approach to the asymp- 
totic value for full QCD. Using the approximate 
value a-^ = 2GeV, they get toq = 1050(170)MeV 
for the quenched theory, and niQ = 730(250)MeV 
for the full theory, where 

2 2,2 /f7\ 

= mo + mg (7) 
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and mg is the degenerate flavour octet {tt/ K/rf) 
mass. 

In addition they have studied the mass spec- 
trum on cooled gauge field configurations and 
have confirmed that the rj' gets its mass from 
large instantons. This agrees with the earlier 
work with Wilson fermions. 

5. SUMMARY AND CONCLUSIONS 

Meson spectroscopy is becoming precise enough 
to ask such questions as whether the quenched ap- 
proximation adequately describes mesons made 
of u, d, s quarks and their antiquarks. Baryon 
spectroscopy is approaching this stage. However, 
understanding of the systematic errors associated 
with choices of sources, sinks, and choices of fits, 
is not yet good enough to preclude significant dis- 
agreements between groups. (From Section 2 one 
would conclude that the quenched approximation 
is less capable of reproducing the correct light 
hadron spectrum than was claimed by the GF- 
11 group §1.) 

With Wilson quarks there does not appear to 
be any strange quark mass which gives the correct 
light hadron spectrum in the quenched approx- 
imation, at the lattice spacings currently used. 
Going to cloverleaf quarks does not appear to 
improve the situation, and even using the Lep- 
age improved gauge and quark actions does not 
lead to a significant improvement. All these cases 
yield values of J in the 0.36-0.40 range, whereas 
the experimental value is in the range 0.48-0.50. 
Preliminary results for staggered quarks extrap- 
olated to a = give J = 0.43(1). The splittings 
in the baryon multiplets (Wilson quarks) appear 
to be too small, for any reasonable value of the 
strange quark mass. 

With good statistics it appears possible to ac- 
curately extrapolate to the "physical" u and d 
quark masses. Such mass extrapolations, and the 
interpolations used to include the effects of the 
s quark, are helped by the fact that, at least for 
light quark physics, the assumption that hadron 
masses depend on only the sum of valence quark 
masses appears valid. 

For full QCD, cloverleaf quarks are an improve- 
ment over Wilson quarks. The value of J is better 



for cloverleaf than Wilson quarks, and for both 
Wilson and cloverleaf quarks it is better than the 
quenched values. 

We have seen the first attempts at implement- 
ing the Lepage scheme for improved actions. For 
more serious implementations this will require go- 
ing to an 0{a'^) fermion action (to match the 
gauge action). This scheme shows great promise 
as a calculational tool. 

For the chiral limit of quenched lattice QCD 
with staggered quarks, it is difficult to sort out 
chiral logarithms from finite size effects and finite 
lattice spacing effects. If we have not seen chi- 
ral logarithms, the obvious question is where are 
they, and what would it take to see them. Part 
of the difficulty in observing chiral logarithms is 
the flavour symmetry breaking, which gives the 
non-Goldstone pions and the disconnected part 
of the 77', larger masses than the Goldstone pion. 
It would therefore be helpful to study quenched 
QCD with Wilson/clover leaf quarks, which have 
no flavour symmetry breaking, down to such light 
pion masses. In conclusion, although the diffi- 
culty in observing quenched chiral logarithms is 
bad news for theory, it means that their effects are 
small, which is good news for phenomenology. 

Although we are still quite a way from doing 
precision calculations of the light hadron spec- 
trum with dynamical quarks, progress has been 
made, particularly within the quenched approxi- 
mation. Improved actions probably are our best 
chance for full QCD. More work appears neces- 
sary to understand systematics. 

Note added in proof: Since the preparation of 
this manuscript we have received a revised draft 
of the LANL preprint In this, a more ex- 

tensive analysis of the baryon octet mass split- 
tings reveals that these splittings are not linear 
in the quark mass splittings, as had been pre- 
viously assumed. When this is properly taken 
into account, the baryon octet splittings are in 
reasonable agreement with experiment. However, 
for the baryon decuplet, the linear formula works 
well, and the disagreement with experiment re- 
mains. 
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